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Gauge fixing to a smooth (Coulomb or Lan- 
dau) gauge on the lattice is essential for the 
extraction of non-perturbative, gauge-dependent 
observables. It also allows a direct comparison 
of non-perturbative results with continuum, per- 
turbative ones. Finally it is emerging as a pow- 
erful approach to measuring gauge-invariant ob- 
servables: like blocking or cooling, it removes 
the noise of ultra-violet fluctuations; but unlike 
these, it preserves the Yang-Mills action. For 
all these reasons, one is interested in generat- 
ing smooth, continuum-like gauge fields on the 
lattice. However the faithful representation of 
continuum gauge fields is a subtle issue, com- 
plicated both by the compact formulation and 
the lattice representation of the underlying space- 
time which allows violations of cherished geo- 
metric identities. For instance, the compact for- 
mulation allows deviations of the Bianchi iden- 
tity dF = d'^A = which introduce monopoles 
and cause lattice QED to confine [1]. Further- 
more, the lattice representation of a fiber bundle 
with the gauge group attached to sites, appar- 
ently has quite a different de Rham cohomology 
than the continuum version. This is responsible 
for the 'lattice Gribov problem' in which gauge 



fixing algorithms find a huge number of config- 
urations which satisfy the gauge condition, in- 
troducing Aharonov-Bohm type objects such as 
Dirac sheets, strings, and vortices. How then do 
we sample over the set of copies, of continuum or 
lattice origin, and what is the proper weight ? 

1. Preliminaries 

To exhibit the deviations of the lattice gauge 
field topology from the continuum which occur 
in gauge fixing, we use the concise geometric no- 
tation of differential forms and here review some 
definitions [3]. A fundamentally useful relation 
is Hodge Decomposition which states that every 
r-form ( in an n-dimensional manifold M) is 
uniquely decomposed into 

where d is the exterior derivative, 6 = (— )"''+"+i * 
d* is its dual via the Hodge star *, and is 
a harmonic r— form that satisfies dh = 6h = 0. 
For a vector field in 3 dimensions, (1) reads in 
coordinates: 

Ai = eijkVjx'' + Vi<l> + hi (2) 
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A more technical point also due to Hodge is the 
relationship 

dim{ = dimH''(M) = K (3) 

between the harmonic forms and the de 

Rham cohomology groups (and hence the 
Betti numbers hr). Recall that classifies, 
roughly, the non-contractible r-spheres in M and 
hence its (global) topology, while are the 

r-dimensional zero eigenvectors of the Laplacian, 
a (local) differential operator. 

1.1. Existence of copies in the continuum 

For a given gauge group G and (compact, Eu- 
clidean) spacetime manifold M , we know when 
to expect generic Gribov copies [4,5]. Consider 
M = X 5"= X ... X 5'"^ a direct product of 
spheres of various dimensions, which includes the 
usual choice of a 4-torus. Then a sufficient condi- 
tion for the existence of Gribov copies is that an 
integer m can be found in {1, 2, ^n^} such 
that the homotopy group i^miG) is non-trivial. 
This theorem arises because gauge transforma- 
tions map the de Rham cohomology structure of 
M (which is rather simple for spheres) onto the 
gauge group. 

For example, all SU(N) groups have 7r3(G) = 
Z, and thus Gribov copies will be found on any 
sphere or torus of c? > 3, including the original 
case studied by Gribov [6]: Coulomb gauge on 
the 3-sphere. 

For G = U (1), 7ri(G) = Z and thus one expects 
copies on any torus (5'"'^ x anything), but not on a 
sphere. The difference is clear: a closed Polyakov 
loop around a torus may have non-trivial wind- 
ing; on a sphere such a loop is contractible and is 
unwound in the process of gauge-fixing. 

1.2. Lattice artifacts 

On the lattice additional gauge copies exist be- 
sides the continuum copies considered above. In 
addition to the large scale (continuum) topology 
of the manifold M , there is a local structure due 
to the elementary lattice hypercubes (plaquette, 
cube, or 4-d hypercube), which allows for non- 
contractible gauge transformations, and gives rise 
to gauge ambiguities. 

These ambiguities are localized lattice artifacts, 
like XY- vortices, Dirac strings, etc... They ap- 



pear because the gauge transformation between 
two lattice 'Gribov' copies generically maps the 
boundary of some hypercubes onto G with non- 
trivial winding [2], and is also a zero mode of the 
lattice Laplacian making it an element of {/}('')}. 
Thus when we construct a lattice fiber bundle by 
attaching a group manifold to each site, the un- 
derlying cohomology of the spacetime appears as 
one with holes at every hypercube, ie. dim if 
= dim {/}'-''-'} diverges in the continuum limit for 
certain r. More accurately it is the homotopy 
groups TTriM) of the fiber bundle which classify 
gauge transformations that are affected. 

Since they do not induce any net global topol- 
ogy (otherwise they would belong to continuum 
copies), these lattice artifacts must come by pairs 
with opposite winding, ie. vortex pairs. The 
gauge transformation which creates a vortex pair 
creates the analog of an Aharonov-Bohm fiux 
tube, whose gauge field is harmonic. However, 
unlike a continuum fiux tube, such a gauge trans- 
formation on the lattice can be continuously de- 
formed into the identity. 

Including these spurious copies in the partition 
function is then overcounting; moreover observ- 
ables may take different values on such copies. 
Thus it is crucial to control the occurrence of vor- 
tex pairs and any bias they may cause. 

2. Numerical gauge fixing by the standard 
(local) method 

The standard approach, in essence, is to ignore 
the Gribov problem: gauge-fixing is performed by 
a local iterative procedure which fails to eliminate 
lattice artifacts, and each gauge-fixed configura- 
tion is given the same weight. The hope is that 
this approach introduces no bias, and converges 
to the correct result in the continuum limit. 

2.1. U(l) on S^: a study of lattice artifacts 

We have examined this expectation for the 
standard choice of lattice Landau gauge 

1 ^ 

f(G,{U}) = -J2ReTrU,'' (4) 

iV — max 

8 = 1 

and for the standard local iterative methods 
of maximizing / (relaxation and over-relaxation 



Figure 1. Vortex pair distribution for cubes of 
size 6 to 12. 

of point singularities, compared to its continuum 
gauge-fixed counterpart. 

Notice that this conclusion could not be 
reached if one would simply monitor the trace 
defect A. As Table I indicates, < A > seems 
to slowly decrease to zero as /3 ^ cxd, mislead- 
ing one toward optimism. In fact, f(G,{U}) is 
the Hamiltonian of a spin-glass with couplings U , 
and the situation we encounter is typical of the 
search for a ground-state, eg. by simulated an- 
nealing: as the system size increases, the energy 
per site approaches that of the true ground state, 
even as the distance in phase space between the 
2 states diverges. 



cube size 


6 


8 


10 


12 


/? 


5 


8.89 


13.89 


20 


<^ vortex pairs> 


0.65 


0.99 


1.09 


1.33 


<trace defect> 


0.017 


0.017 


0.014 


0.013 



Table 1 



Evolution with /3 of the average properties of 
gauge copies. 

To understand better the behavior of A, we 
looked at its correlation with the number of vor- 
tex pairs. Within our statistics, A grew in pro- 
portion to the number of pairs, confirming our 
expectation that vortices form a dilute, non- 
interacting gas. Further evidence was provided 
by the distribution of the mutual angle, on the 
sphere, between a vortex and an anti-vortex: it 
was uniform, except for a sharp drop to zero at 
small angles, corresponding to ~ 2 lattice spac- 
ings. This uniform distribution refiects the failure 
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of the local maximization used, to bring together 
and annihilate vortex and anti-vortex unless they 
are next to each other on the lattice. 

The trace defect per vortex pair can be approx- 
imately calculated by considering the prototypi- 
cal lattice copy on the sphere: a vortex at the 
North pole and anti-vortex at the South pole, as 
if created by a Dirac string piercing the sphere 
along the polar axis. The harmonic gauge field 
produced by the string is A = dcfj/rsinO, and the 
trace defect is „,."'" , [„„ A A *A. Taking for the 
number of links SirR'^/a'^, and integrating over 
the whole surface of the sphere except for 2 cir- 
cular caps of radius a, one gets 



A 



1 



(5) 



which goes to zero with a. In actuality the vor- 
tices are not systematically at opposite poles, and 
the integration is really over a discrete cube, so 
that the coefficients in (5) will be changed. Still 
a formula of the type (5) describes our data well, 
as Fig. 2 shows. 

0.035 r^, , . , , . , 
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0.005 



f(x) = -a * log{4/x**2)/x**2 
a = 0.421 (3) 



8 10 12 14 16 18 



Figure 2. Trace defect per vortex pair, as a func- 
tion of the lattice spacing. 

This simple toy model shows the fatal disease of 
the standard local approach to gauge fixing: as 
the lattice spacing decreases, gauge-fixed 
configurations have less and less to do with 
their continuum counterpart. We do not ex- 
pect this conclusion to change with the gauge 
group or the lattice manifold. 

2.2. Measurable effects of lattice artifacts 

Several attempts have been made at detecting 
observable effects of spurious gauge copies [8-10]. 



This is a very difficult task, since one does not 
know how to systematically identify and elimi- 
nate these copies. At best, one can establish a 
correlation between the value of some physical ob- 
servable and the "quality" of gauge fixing [11,12]. 
It also seems to depend on the sensitivity of the 
observable to topology. 

We decided to look at this issue again on a toy 
model, U(l) on a 2-dimensional torus. Specifi- 
cally we tried to extract the correlation length ^ of 
2 Polyakov loops. In this model each link U is de- 
fined by an angle ^{x) = jXogU ^{x) E [— 7r,7r]. 
We measured this correlation in 3 different ways: 

(1) < e'(®^(*)-®^(°)) > This is the usual 

gauge-invariant procedure, where Qp{t) = 

'^x=i^^i^'^)- ^® repeated the measure- 
ment for 3 values of [3, from 16 to 256, on 
16^ to 64^ tori, keeping the physical size 
constant. The lattice model is exactly solv- 
able in this case [13], providing a check of 
our Monte Carlo program. 

(2) < {Qp{t) - ep(0))2 >. This is the correla- 

tion one would measure after fixing to the 
nearest maximum of / in Coulomb gauge, 
with ex{x,t) = l/L ep(t). 

(3) < {ep{t) - Qp{0)f > where 6 is the pro- 

jection of in [— TT, tt]. This corresponds to 
fixing to Coulomb gauge and selecting the 
global maximum of /. 

One can see in Fig. 3 that (3) is a noisy but 
correct way to measure the correlation length ^, 
but that (2) is completely wrong. In addition, 
there is no trend for (2) to approach the correct 
result as a decreases. Thus we have exhibited one 
situation where lattice copies completely spoil a 
physical measurement, and appear to continue to 
do so in the continuum limit. 

3. Gauge fixing by Hodge decomposition: 
a global approach 

3.1. Abelian gauge fields 

Since we have shown above that the lattice Gri- 
bov problem is essentially due to lattice artifact 
harmonic forms, Hodge decomposition provides 



5 



0.5 



0.2 



0.1 



0.05 



0.02 



< A,(0) A,(t) > 



10 20 30 

Figure 3. Gauge-fixed correlations on tlie torus, 
witli (daslied line) and without (dotted line) spu- 
rious gauge copies. The solid line is the gauge- 
invariant analytic result. 

us with a very powerful tool for separating the 
continuum gauge fixed fields from the lattice de- 
fects. While the exterior algebra on a lattice [14] 
is defined for non-compact complex fields, we ex- 
tend it to an abelian lattice fiber bundle by using 
link angles 9fj,{x) = Afi(x). From (1) we have 
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Ax (6) 
= Act,. (7) 
) = 6d is the lattice Laplacian 

(^Aj^'j^ — ^^^^ ^x'x-\-ii ^^x'x ~1~ ^x'x — iJ,)- 

The Landau gauge condition is 6 A = 0, which 
when translated into links is 



dA 
6A 
Here A 



6A 



1 'ix. 



(8) 



We take this as the definition of our Landau gauge 
condition for abelian fields, and will address non- 
abelian fields shortly. 

There are 2 possible methods to fix to Landau 
gauge then. The first one is to compute 6A(x) = 
^^[j4^(a;) — Afi(x — fi)]. From this compute (t)(x) 
in the Hodge decomposition of A by inverting the 
Laplacian in (7) (via the Fast Fourier Transform). 
From (j), make the gauge transformation 

A' = A - d<j, = 6x + h. (9) 

A' now satisfies 6A' = identically, but is not 
unique since h may be any one of the lattice ar- 
tifact harmonic forms. 



We could proceed and eliminate the harmonic 
part h by finding the gauge transformation cor- 
responding to the topological object described by 
h (eg. vortex pairs in 2-D), however this is cum- 
bersome. It is far simpler to just set 



A^A' = 6x = 6(A-^F). 



(10) 



We display this method in Figs. 4-6 (L=16). The 
initial field Ui(x) = -I sin(27ra;/_L) cos(27rt//_L), 
U'zix) = 1, was transformed with a random gauge 
transformation (Fig. 4). The gauge fixed config- 
uration is shown in Fig. 5. For comparison, we 
show also the same field gauge fixed by the local 
method of section 3 in Fig. 6. 
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Figure 4. Initial gauge field configuration. 
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Figure 5. Landau gauge fixed field configuration: 
Hodge method A' = (5(A-i_F). 

Notice that we have some freedom in comput- 
ing F in (10) which is very interesting. If we 
use the compact definition of F = j log[/piq, the 
gauge fixed field A' will preserve the DeGrand- 
Toussaint monopoles which violate dF = 0. Us- 
ing the non-compact definition of i^^^ = A^(x + 
fi) — A^(x) — Afi(x + ly) + Afi(x), will set A' to 
the equivalent non-compact gauge field configu- 
ration, without lattice monopoles. Thus we have 
a way of easily converting a field from compact 
form to non-compact form, preserving the action, 
and gauge fixing at the same time. 
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Figure 6. Gauge fixed field configuration: Stan- 
dard maximization method of section 3 

3.2. Non-abelian gauge fields 

The above gauge fixing method works perfectly 
for abelian fields since the discrete exterior alge- 
bra lifts (up to compactness issues) nicely to the 
abelian fiber bundle. However, the situation is 
rather complicated for non-abelian fields. We still 
have Hodge decomposition (1) for Lie algebra val- 

(r) 

ued r— forms Aa r", however dA is not a gauge 
invariant quantity, since it is the covariant exte- 
rior derivative of A, DA = F which is the field 
strength. Furthermore, nilpotence is lost for D 
on Lie algebra valued forms: _D^$ = F A^. Thus 
the program above must be modified. 

We found interestingly, that if we apply the 
first method outlined for the abelian case, ie. ex- 
tract (j>" = A~^(6A") and make the gauge trans- 
formation U G\x)Uij,{x)G{x + fi), where a is 
the SU (2) algebra index and G = exp(— iec/)) that, 
for < e < 1.3, successive gauge transformations 
iteratively converge to 6 A" = 0. 

There are a variety of geometrical extensions of 
this method which we are exploring. We present 
next a simple, smooth, and unique gauge for non- 
abelian fields. 

3.3. The r gauge 

Since we have a unique gauge fixing method for 
abelian fields we propose the following algorithm 
for gauge fixing non-abelian fields. 

First, choose an operator 0(x) which trans- 
forms equivariantly (0'^{x) = G\x)0{x)G{x)) 
and can thus be rotated at each lattice site to a 
fixed direction in V space by a gauge transforma- 
tion. This leaves only an abelian gauge symmetry 
corresponding to transformations which commute 
with the chosen vector in r space. Then, follow- 
ing the above method based on Hodge decompo- 



sition, gauge fix this remaining abelian freedom. 

For an SU(2) example, take O to be the sum 
of plaquettes at each site, projected onto SU(2), 
and diagonalize it. We are left free to make gauge 
transformations G(x) = exTp[ia(x)Ts\ which we 
can use to set, say, 6A^ = 0. Thus this gauge has 
a continuum limit of 

^i^,V = 0, 5,4 = (11) 

This gauge could be a simple and useful choice, 
especially for the study of monopole-induced con- 
finement since it is very much like t'Hooft's orig- 
inal proposal. It is economical, since the global 
step can be performed by FFT. It is unique, up to 
known continuum harmonic forms. Furthermore 
it has a convenient perturbative interpretation, 
unlike earlier global proposals [15]. 
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